METHOD OF AND APPARATUS FOR 
SEARCHING CORRESPONDING POINTS BETWEEN IMAGES, AND 
COMPUTER PROGRAM 

5 FIELD OF THE INVENTION 

The present invention relates to a technology of 
searching corresponding points for establishing 
correspondence between an input image and a reference image 
which is an obj ect for comparison with the input image . More 
10 particularly, this invention relates to a technology capable 
of rapidly obtaining a stable result without resulting in 
obtaining a local solution ( local minimum) when establishing 
a correspondence between the input image and the reference 
image . 

15 

BACKGROUND OF THE INVENTION 

Conventionally, a method of establishing 

correspondence between images by searching corresponding 

points between an input image and a reference image when 
20 verifying the input image with a reference image beforehand 

registered. The input may be input through an image input 

apparatus such as a scanner. 

For example, in "Computational vision and 

regularization theory" by T.Poggio, V.Torre and C.Koch, in 
25 NATURE Vol.317, pp. 314-319, 1985 ("conventional art 1"), 
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discloses a technique using the standard regularization 
theory for minimizing a set energy function in a calculus 
of variations. In the standard regularization theory 
adopted in the conventional art 1, corresponding points 
5 between images are calculated at minimizing energy by 
employing a repetition calculation merely using a local 
information. As a result, parallel distributed process can 
be performed, and in addition, information processing such 
as that performed in a human brain may be realized. 

10 "Monotonic and continuous two-dimensional warping 

method based on a Dynamic Programming " by Seiichi UCHIDA 
and Hiroaki SAKOE, in THE TRANSACTIONS OF THE INSTITUTE OF 
ELECTRONICS, INFORMATION AND COMMUNICATION ENGINEERS D-II, 
Vol. J81-D-II no . 6, pp . 1251-1258 , June 1998 ("conventional 

15 art 2"), discloses a technique of a two-dimensional 
DP-warping method for efficiently searching an optimum 
solution with DP. According to the conventional art 2, the 
optimization problem can be solved efficiently. 

However, the conventional art 1 has a problem . Because 

20 the corresponding points are calculated by minimizing energy 
with repetition calculation me rely using a local information, 
the solution greatly depends on the initial value. In 
addition, because it is easy to converge to a local solution, 
it is difficult to obtain an optimum correspondence. 

25 Further, the conventional art 2 also has a problem. 



Although an optimum solution can be efficiently searched 
with DP, amount of calculation becomes vast. Concretely, 
calculation time having exponent-order for an image size 
is required in order to obtain the optimum solution. 
5 In view thereof, an important problem occurs how an 

apparatus is realized which searches corresponding points 
of images and which can rapidly obtain a result of a stable 
correspondence without resulting in any local solution when 
making an image correspond to the other image, i.e. applying 
10 correspondence relationship to between images. 

SUMMARY OF THE INVENTION 

It is an object of this invention to provide a method 
and an apparatus capable of rapidly obtaining a stable result 

15 without resulting in obtaining a local solution (local 
minimum) when establishing a correspondence between the 
input image and the reference image. It is another object 
of this invention to provide a computer program containing 
instructions which when executed on a computer realizes the 

20 method according to the present invention. 

According to the present invention, a plurality of 
similarity degree images each having a similarity degree 
between the input image and the reference image as a pixel 
value are produced. Then, corresponding points between the 

25 input image and the reference image are detected based on 
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the detected plurality of similarity degree images 
Other objects and features of this invention 
become apparent from the following description 
reference to the accompanying drawings. 

5 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig.l shows a functional block diagram showing a 
construction of a corresponding point searching apparatus 
relating to the embodiment of the present invention. 
10 Fig.2A is an example of a reference image and Fig. 

2B is an example of an input image used in the embodiment 
of the present invention. 

Fig. 3 shows an example of a reference partial image 
(i.e. a portion of the reference image). 
15 Fig. 4 shows examples of input partial images (i.e. 

portions of the input image). 

Fig. 5 explains how a similarity degree image is 
prepared . 

Fig. 6 shows a similarity degree image corresponding 
20 to each of the input partial image shown in Fig. 4. 

Fig. 7 explains an accumulation-addition procedure. 

Fig. 8 explains the accumulation-addition process in 
a j-direction. 

Fig. 9 explains the accumulation-addition process in 
25 the -j-direction (opposite of the j-direction). 



Fig. 10 explains the accumulation-addition process in 
a i-direction. 

Fig. 11 explains the accumulation-addition process in 
the -i-direction (opposite of the i-direction) . 
5 Fig. 12 shows similarity degree images after 

accumulative calculation. 

Fig. 13 shows positions of maximum value in the 
similarity degree images shown in Fig. 12. 

Fig.l4A and Fig. 14B show examples of distortion 
10 (deformation) of the reconstructed input image based on the 
corresponding points. 

Fig. 15 is a flow chart of a process performed by the 
corresponding point searching apparatus shown in Fig.l. 

Fig. 16 is a flow chart of an example of the 
15 accumulation-addition process. 

Fig. 17 is a flow chart of an example of the 
accumulation-addition process when a maximum value filter 
is used. 

Fig. 18 is a flow chart of process performed by the 
20 maximum value filter. 

Fig. 19 explains how a character may be cut off from 
the input image consisting of a many characters. 

Fig. 20 shows an example in which the present invention 
is applied to a recognition of a Chinese character. 
25 Fig. 21 shows an example in which the present invention 



is applied to signature authentication. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Preferred embodiment of a method of and an apparatus 
5 for searching corresponding points between images to search 
corresponding points for making an input image correspondent 
to a reference image of object of comparison with the input 
image, and a computer program for realizing the method 
according to the present invention on a computer will be 

10 explained below with reference to the accompanying drawings . 

A functional block diagram of the apparatus for 
searching corresponding points according to the present 
embodiment is shown in Fig. 1. This apparatus 1 mainly 
comprises (1) a similarity degree image production section, 

15 and (2) a corresponding points determination section . The 
similarity degree image production section produces a 
similarity degree image for indicating a similarity degree 
between an input partial image obtained by dividing the input 
image and the reference partial image obtained by dividing 

20 the reference image. The corresponding points 

determination section determines optimum corresponding 
points based on the result obtained by accumulatively adding 
a plurality of similarity degree images. 

Concretely, in the similarity degree image production 

25 section, e.g. a correlation value is sequentially obtained, 

6 



moving a reference partial image formed by 7x7 pixel size 
in an input partial image formed by 21 x 21 pixel size and 
a similarity degree image having this correlation value as 
a pixel value is produced. As a result, this similarity 
5 degree image will indicate to what extend a central point 
on the reference partial image is similar to a point on an 
input partial image. 

In the corresponding points determination section, 
each of similarity degree images is renewed, recursively 

10 repeating accumulation-addition for a pixel value of each 
of the similarity degree images in a vertical direction and 
in a horizontal direction to produce a similarity degree 
image and obtains a position of pixel having a maximum image 
value of images forming each similarity degree image is 

15 obtained, so that this position is determined as a 
corresponding point on the input partial image. 

Thus, according to the apparatus 1 for searching 
corresponding points, it becomes harder to result in 
obtaining a local minimum because a wider optimization can 

20 be performed as compared to the conventional technique. 
Correspondence between images becomes highly accurate 
because an optimizing minute adjustment repeating 
accumulation-addition is employed. And a hardware becomes 
easy to fabricate because a parallel distributed process 

25 becomes possible. 



Construction of the apparatus 1 for searching 
corresponding points shown in Fig.l will be explained in 
more detail. It will be considered here that images are 
input in the order of reference image and input image. As 
5 shown in Fig . 1 , the apparatus 1 comprises image input section 
10, division process section 11, reference partial image 
temporal storage section 12, similarity degree image 
production section 13, accumulation-addition processing 
section 14, and corresponding points determination section 
10 15 . 

The image input section 10 receives a reference image 
Io(i, j) and an input image Ii(i, j). Here, O^i^I-1 and 
0^= j J-l) , and I is a number of pixels in vertical direction 
and J is a number of pixels in horizontal direction. 

15 Concretely, the image input section 1 0 may be an image scanner 
which optically read an object to obtain an image. On the 
other hand, image input section 10 may be an interface section 
for obtaining an image from a network, or it may be a reading 
section for reading an image from a memory and so on. The 

20 input image is defined as an image of object to be searched 
for corresponding points and it may accompany with distortion 
and deformation. On the contrary thereto, the reference 
image is defined as an image to be compared with the input 
image and preferably, it may not have distortion or the like . 

25 Fig.2A shows the reference image 21. Fig.2B shows 



the input image 22. Both, the reference image 21 and the 
input image 22, are images of English alphabet "a". Both 
the images are of size 32 x 32 pixels. 

The division process section 11 is for dividing an 
5 input image and a reference image input from the input image 
portion 10 into an input partial image and a reference partial 
image, respectively. Here, the procedure of division of 
the input image and that of the reference image are mutually 
different . 

10 When dividing the reference image, a reference partial 

image is produced, whose center lies at a center of sampling 
point (p m , q n ) (where O^m^M-1, O^n^N-1) obtained by 
sampling M numbers of points in a vertical direction and 
N numbers of points in a horizontal direction on the reference 

15 image. 

Fig. 3 shows an example of result of division of the 
reference image 21 shown in Fig.2A. As shown in Fig. 3, the 
reference image 21 formed by 32 x 32 pixels is divided into 
twenty five reference partial images formed by 7 x 7 pixels . 

20 Concretely, the relationship of p m = round ( I/M ) and q n 
= round ( J/N ) is satisfied. The term of "round ( ) " is meant 
by "round to the nearest whole number". When dividing an 
input image, being different from dividing the reference 
image, the input image is divided so as to have overlapped 

25 data in which parts of each of the input partial images are 



mutually overlapped. 

Fig. 4 is for showing an example of result of division 
of the input image 22 shown in Fig.2B. As shown in Fig. 4, 
the reference image 22 formed by 32 x 32 pixels is divided 
5 into twenty five reference partial images formed by 21 x 
21 pixels. 

The reference partial image temporal storage section 
12 is for temporally storing each of the reference partial 
images divided by the division process section 11 in which 

10 a corresponding reference partial image is picked up when 
the similarity degree image production section 13 produces 
a similarity degree image. 

The similarity degree image production section 13 is 
for calculating a similarity degree, considering over 

15 deformation between the input partial image and the re f erence 
partial image to produce a similarity degree image Cmn (u, v) 
(where O^m^M-1, O^n^N-1, O^u^U-1 and O^v^V-1) having 
the similarity degree as a pixel value. Here, U and V 
respectively corresponds a vertical size and a horizontal 

20 size . A normalized correlation coefficient ( a f g / ( a f a g ) ) 
can be used as this similarity degree. 

Fig .5 is for explaining about conception for producing 
a similarity degree image from the reference partial image 
and the input partial image. Here, the input partial image 

25 and the reference partial image can be used which are 



positioned at a position of two-rows and two-columns in Fig . 3 
and Fig. 4. 

When obtaining a similarity degree between the input 
partial image 51 and the reference partial image 52 shown 
5 therein is obtained, a normalized correlation coefficient 
is calculated, making a center pixel of the reference partial 
image 52 correspond to a pixel at an upper-left portion of 
the input partial image 51. The calculation result is 
defined as apixel value of the pixel at the upper- left portion 

10 of the similarity degree image 53. Thereafter, the 
reference partial image 52 is displaced rightward and the 
process is performed in a same way. The above-mentioned 
process is performed for all the pixels of the input partial 
image 51, displacing the reference partial image 52, 

15 resulting in obtaining the similarity degree image 53 . 

If such a process for producing a similarity degree 
image is performed for each of input partial images, then 
a plurality of similarity degree images can be obtained as 
shown in Fig. 6. Here, if each of pixel values of all the 

20 pixels in the reference partial images is constant, then 
the denominator of a normalized correlation coefficient 
becomes zero. As a result, the pixel value of the similarity 
degree image in this case also becomes zero. 

The accumulation-addition processing section 14 is 

25 for recursively accumulation-adding each of similarity 



images in order from in the j -direction, the -j -direction, 
the i-direction and the -i-direction . Concretely, when 
accumulation-adding in the j -direction for similarity 
degree images of n=l to (N-l) ; 
5 Cmn(u,v)=Cmn(u,v)+ aMAX (Cmn_i (p, q) ) 

which is sequentially and recursively calculated. Here, 
the term of "Max ( )" is defined as a maximum value and a 
is a constant value (where O^u^U-1, O^v^V-1, O^m^M-1, 
O^n^N-1, u-l^p^u+1 and v-l^q^v+1) . 

10 Namely, when performing accumulation-addition in the 

j -direction, a maximum value is obtained in 3x3 pixels whose 
center is C mn-l(p,q) as shown in Fig. 8. Thereafter, this 
maximum pixel value is multiples by a. The multiply value 
is added to the pixel value of Cmn(p,q) . Such a process is 

15 recursively repeated. 

Further, for a similarity degree image of n=(N-2) to 
0 when performing accumulatively addition in - j -direction, 
following equation is sequentially and recursively 
calculated . 

20 Cmn(u,v) = Cmn(u,v) + aMax (Cmn+i (p, q) ) 

where O^v^V-1, and Namely, 

when performing accumulation-addition in the - j -direction, 
as shown in Fig . 9 , a maximum pixel value in 3x3 pixels whose 
center lies at Cmn+i ( p, q ) is obtained. This maximum pixel 

25 value is multiplied by a. Then a process is recursively 



repeated for adding the multiply value to the pixel value 
o f Cmn ( p , q ) . 

Furthermore, for similarity degree images of m=l to 
(M-l) when performing accumulation-addition in the 
5 i-direction, following equation is sequentially and 
recursively calculated; 

Cmn(u,v) = Cmn(u,v) + aMax ( Cm.- in ( p , q) ) 
where Namely, 
when performing accumulation-addition in the i-direction, 
10 as shown in Fig . 1 0 , a maximum pixel value in 3x3 pixels whose 
center lies at Cm_in ( p, q ) is obtained. This maximum pixel 
value is multiplied by a. Then a process is recursively 
repeated for adding the multiply value to the pixel value 
o f Cmn ( p , q ) . 

15 Furthermore, for a similarity degree image of m= (M-2 ) 

to 0 when performing accumulation-addition in the 
-i-direction, following equation is sequentially and 
recursively calculated. 

Cmn(u,v) = Cmn(u,v) + aMax(Cm+in (p,q)) 

20 where O^u^U-1, O^v^V-1, O^m^M-2, and O^n^N-1 . Namely, 
when performing accumulation-addition in -i-direction, as 
shown in Fig. 11, a maximum pixel value in 3x3 pixels whose 
center lies at Cm+in ( p, q ) is obtained. This maximum pixel 
value is multiplied by a. Then a process is recursively 

25 repeated for adding the multiply value to the pixel value 



o f Cmn ( p , q ) . 

The corresponding points determination section 15 is 
for determining corresponding points based on a similarity- 
degree image after accumulation-addition in cooperation 
5 with the accumulation-addition processing section 14. 
Concretely, a position of a maximum value in each of 
similarity images is performed when the 

accumulation-addition is completed in each of the directions 
Thereafter, if variance between the position thereof and 

10 the position of the maximum value at the time of the previous 
accumulation addition does not lie within a predetermined 
area, then the accumulation addition process is feedbacked 
and repeated. The repetition process is completed at the 
time when it lies within the predetermined area. Then the 

15 position of the maximum value of each of similarity degree 
images is decided as a corresponding point. 

For example, when obtaining each of the similarity 
degree images as shown in Fig. 12 by repetition of 
accumulation addition process, an amount of variance of 

20 position of the maximum value of each of the similarity degree 
images is examined. When the amount thereof lies to exist 
within a predetermined area, then a position of the maximum 
value of each of the similarity degree images is obtained 
as shown in Fig. 13. This obtained points are defined as a 

25 corresponding points. Here, distortion of the input image 



based on the corresponding points will be shown in Fig. 14. 

Next, a process performed by the apparatus 1 for 
searching corresponding points shown in Fig.l will be 
explained with reference to Fig. 2 through Fig. 14. 
5 Fig. 15 shows a flow chart for indicating a process 

procedure of the apparatus 1 for searching corresponding 
points shown in Fig . 1 . As shown in Fig. 15, if the apparatus 
1 for searching corresponding points obtains a reference 
image and an input image, then the reference image is divided 

10 into reference partial images (Step 101) . For example, the 
reference image 21 formed by 32 x 32 pixels shown in Fig.2A 
is divided into twenty five reference partial images formed 
by 7 x 7 pixels shown in Fig. 3 and the divided ones are stored. 
This reference partial images may be pre-stored by each of 

15 images to be referred. Not a reference partial image but 
an only reference image is pre-stored independently, the 
stored reference image may be divided into reference partial 
images at the time of below-mentioned calculation of 
similarity between an input image and a reference image. 

20 Thereafter the input image is divided into input 

partial images (Step 102) . For example, the input image 
22 formed by 32 x 32 pixels shown in Fig.2B is divided into 
twenty five input partial images formed by 21 x 21 pixels 
shown in Fig. 4. 

25 Then, a similarity degree image is produced between 



the reference partial images and the input partial images 
(Step 103) . For example, as shown in Fig. 5, the similarity 
degree image 53 is produced in which similarity degrees of 
the input partial image 51 and the reference partial image 
5 52 are defined as a pixel value produced by each of the input 
partial images, so that a plurality of similarity degree 
images as shown in Fig. 6 are produced. 

Next, the j -direction accumulation addition process 
(Step 104) , the -j-direction accumulation addition process 

10 (Step 105), the i-direction accumulation addition process 
(Step 106), and the -i-direction accumulation addition 
process (Step 107) are performed in which a position of a 
maximum value of the similarity degree image is detected 
( Step 108 ) . Concretely, addition processes shown in Figs . 8 

15 to 11 are recursively repeated as shown in Fig. 7, the 
similarity degree image as shown in Fig. 12 is produced and 
a position of the maximum value of each of the similarity 
degree images is detected as shown in Fig. 13. 

Next, itis examined whether or not variance of position 

20 of the maximum value lies within a constant value (Step 109). 
If not so (No at Step 109) , then the procedure is advanced 
to Step 104 in which a same way is repeated. On the other 
hand, if so (Yes at Step 109) , the procedure is completed, 
defining this position as a corresponding point. 

25 Next, an explanation will be given about the 



j-direction accumulation-addition process procedure using 
the accumulation-addition processing section 14 shown in 
Fig.l. Fig. 16 shows a flow chart for indicating the 
j-direction accumulation addition process procedure using 
5 the accumulation-addition processing section 14 shown in 
Fig . 1 . Here , such an accumulation addition process can also 
be applied to the -j-direction, the i-direction, and the 
-i-direction ones. 

As shown in Fig. 16, at first, variables m, u and v 
10 are set as zero and initialization is performed for setting 
the variable n as one {Steps 201 to 204) . Here, the variable 
m is used as an index in the i-direction. The variable n 
isusedasan index in the j -direction index and the variables , 
variables u and v are relatively ones in the i-direction 
15 and the j-direction for indicating a searching area. 

If this initialization is completed, then following 
calculation is performed (Step 205) . 

Cmn(u,v) = Cmn(u,v) + ocMax (Cmn_i (p,q)) 

Thereafter, the variable v is incremented (Step 206) . 
20 If this variable v is smaller than V (Yes at Step 207) , then 
the procedure is advanced to Step 205 in which the addition 
process is repeated. Namely, an area to be searched is 
displaced in the j-direction. 

On the contrary thereto, if the variable v is more 
25 than V (No at Step 207) , then the variable u is incremented 



(Step 208) . If the variable u is smaller than U (Yes at Step 
209) , then the procedure is advanced to Step 204 in which 
addition process is repeated . Namely, an area to be searched 
is displaced in the i-direction. 
5 If the variable u is more than U (No at Step 209), 

then calculation for a pixel is completed and the procedure 
is advanced to a next pixel. Concretely, the variable n 
is incremented and a target pixel is advanced in the 
j-direction (Step 210) . Thereafter this variable n is 

10 compared with N (Step 211) . If the variable n is smaller 
than N (Yes at Step 211), then the procedure is advanced 
to Step 203 and the addition process is repeated. 

On the contrary thereto, if the variable n is more 
than N (No at Step 211) , then the variable m is incremented 

15 and a target pixel is advanced in the i-direction (Step 212 ) . 
If this variable m is smaller than M (Yes at Step 213) , then 
the procedure is advanced to Step 2 02 and the addition process 
is repeated. If the variable m is more than M (No at Step 
213) , then the process is completed. According to a series 

20 of process as above, a result of accumulation-addition in 
the j-direction for all the pixels in each of similarity 
images can be obtained. 

By the way, a maximum value is calculated, incrementing 
variables u and v in the flow chart of Fig. 16 . As the other 

25 calculation method, a maximum value filter image C mn_i (u, v) 



to which a maximum value filter is applied for Cmn_iis produced 
beforehand and thereafter it is multiplied by a. The 
multiply result may be added to Cmn(u,v). 

Here, a process procedure when using such a maximum 
5 value filter will be hereinafter explained using Fig. 17 
and Fig. 18. Fig. 17 shows a flow chart for indicating the 
j-direction accumulation-addition process procedure using 
a maximum value filter with the accumulation-addition 
processing section 14 shown in Fig.l. Fig. 18 shows a flow 
10 chart for indicating a process procedure of the maximum value 
filter . 

As shown in Fig. 17, the variation m is initialized 
to be "0" (Step 301) and the variation n is initialized to 
be "1" (Step 302) . Thereafter the maximum value filter 

15 Cmn_i(u,v) is calculated (Step 303). 

Concretely, as shown in Fig. 18, the variations u and 
v are initialized to be zero (at Steps 401 and 402) . 
Thereafter, the equation of C'mn(p,q) = Max [Cmn(p,q)] is 
operated (Step 403 ) . Here, the condition of u-l^p^u+1 

20 and v-l^q^v+1 is satisfied. The variable v is incremented 
( Step 404 ) . if the variable v is smaller than V (Yes at 
Step 405), the procedure is advanced to Step 403 and 
calculation as above is performed. If the variation v is 
more than V (No at Step 405), then the variation u is 

25 incremented (Step 406) . If the variation u is smaller than 



U (Yes at Step 406) , the procedure is advanced to Step 402 
and process such as above is repeated, so that a maximum 
value is obtained. 

Thus, if a maximum value filter is defined through 
5 the calculation, then the variables u and v are initialized 
to be "0" (Steps 304 and 305) . Thereafter, the following 
equation will be operated (Step 306) : 

Cmn(u,v) = Cmn(u,v) + aMax (C'mn-i (p,q)) 

Thereafter, the variation vis incremented (Step 307 ) . 
10 If the variation v is smaller than V (Yes at Step 308) , then 
the procedure is advanced to Step 306 and addition process 
is repeated. Namely, an area to be searched is displaced 
in the j-direction. 

On the other hand, if the variable v is more than V 
15 (No at Step 308) , then the variable u is incremented (Step 

309) . If the variable u is smaller than U (Yes at Step 310) , 
the procedure is advanced toStep 305 and the addition process 
is repeated. Namely, an area to be searched is displaced 
in the i-direction. 

20 Further, if the variable u is more than U (No at Step 

310) , then calculation for a pixel is completed and the 
procedure is advanced to a next pixel. Concretely, the 
variable n is incremented and the target pixel is moved in 
the j-direction (Step 311) . When comparing this variable 

25 n with N (Step 312), if this variable n is smaller than N 



(Yes at Step 312), then the procedure is advanced to Step 
303 and the addition process is repeated. 

On the contrary thereto, if the variable n is more 
than N (No at Step 312), the variable m is incremented and 
5 the target pixel is moved in the i-direction (Step 313) . 
Thereafter, if the variable m is smaller than M (Yes at Step 
314), then the procedure is advanced to Step 302 and 
accumulation process is repeated. If the variable m is more 
than M (No at Step 314), then the process is completed. 
1° Next, an explanation will be given about a process 

after searching corresponding points . If the corresponding 
points are searched using the corresponding point searching 
apparatus 1 shown in Fig.l, then it will be apparent how 
the input image is distorted, comparing the input image with 
15 the reference image in view of these corresponding points. 

Concretely, as shown in Fig. 19, when de fining an image 
having a numerical sequence of "523" as an information, as 
an input image and defining an image having a numeral of 
"2" as an information, as a reference image, it can be 
20 understood that distortion occurs as shown therein. 
Additionally, when defining an image having a numeral of 
"3" as an information, as a reference image, it can be 
understood that distortion occurs as shown therein. 

Like this, distortion is used based on such 
25 corresponding points. Thereby, each character is 



separated from a continuously written character sequence, 
so that an accurate character recognition becomes possible. 

This point will be further explained hereinafter . For 
example, when selecting an image having a numeral of "0" 
5 as a reference image, a portion of a closed curve between 
the numerals of "5" and "2" in the input image is defined 
as "0", resulting in determining corresponding points, i.e. 
the distortion. The degree of this distortion is greater 
than that of the degree when defining an image having a numeral 

10 of "2" as a reference image, so that it can be decided that 
there is no numeral of "0" in the input image. 

Further, as shown in Fig. 20, when defining an image 
having a Chinese character meant by "bird" as an information, 
as a reference image and an input image, the distortion as 

15 shown therein occurs in the input image. Accordingly, the 
similarity degree is calculated, reforming the distortion, 
i.e. the deformation, so that an accurate similarity degree 
can be grasped resulting in that verification of accuracy 
can be improved. 

20 Furthermore, as shown in Fig . 21 , the present invention 

embodiment can also be applied to signature authentication. 
It is very difficult to authenticate a signature drawn in 
English. However, if the degree of the distortion based 
on corresponding points is obtained using the present 

25 invention, then a registration sign and an authentication 



sign to be authenticated can be accurately authenticated. 
Here, for the sake of explanation, though drawings will be 
omitted, the present invention can be applied to detection 
of an optical flow, 3D-measurement using a stereo image, 
5 and an apparatus for examining an external appearance and 
so on . 

As above-mentioned, according to the present 
embodiment, a reference image and an input image are input 
from an image input section 10. The division process section 

10 11 produces a reference partial image and an input partial 
image. The similarity degree image production section 13 
respectively produces similarity degree images whose 
similarity degree between the reference partial image and 
the input partial image. The accumulation-addition 

15 processing section 14 recursively performs an 
accumulation-addition of a plurality of similarity degrees . 
Based on the result, the corresponding points determination 
section 15 determines corresponding points. Thus, (1) a 
case of resulting in obtaining a local minimum can be reduced 

20 because optimization is widely performed, (2) the 
optimization can be minutely adjusted using a repetition 
process, and (3) a rapid process can be realized because 
a parallel distributed process can be performed. 

Namely, e.g. in the above embodiment, for the sake 

25 of explanation, a case is shown where the size of a reference 



image is as same as that of an input image. The present 
invention is not limited by such an embodiment . For example, 
this invention can also be applied to a case where the size 
of a reference image is different from that of an input image . 
5 Additionally, in the above embodiment, a case is shown 

where each of similarity degree images is recursively 
accumulation-added in the order of j -direction, 
-j -direction, i-direction, and -i-direction . The present 
invention is not limited by such an embodiment, too. For 

10 example, this invention can also be applied to a case where 
each of similarity degree images is accumulation-added in 
an oblique direction or an only single direction. Even a 
case of in the only single direction, similarity degrees 
are added, so that a stable information can widely be 

15 obtained. 

Further additionally, in the above embodiment, a case 
is shown where a normalized correlation coefficient is used 
as a similarity degree . The present invention is not limited 
by such an embodiment. 

20 As explained above, this invention provides an 

apparatus capable of rapidly processing with a repetition 
process by each of blocks. 

Furthermore, this invention provides an apparatus 
capable of performing a process considering over a distortion 

25 (deformation) of an input image. 



Furthermore, this invention provides an apparatus 
capable of obtaining a similarity degree image being robust 
for variance of contrast of an image. 

Furthermore, this invention provides an apparatus 
5 capable of widely using an information and reducing a case 
resulting in obtaining a local minimum. 

Furthermore, this invention provides an apparatus 
capable of performing a process in which information is used 
more widely. 

10 Furthermore, this invention provides an apparatus 

capable performing an optimized minute adjustment by a 

repetition process. 

Furthermore, this invention provides an apparatus 

capable of easily specifying corresponding points. 
15 Furthermore, this invention provides an apparatus 

capable of saving useless accumulation-addition and 

efficiently correcting a deformation. 

Furthermore, when applying correspondence 

relationship to between images, this invention provides a 
20 method capable of rapidly obtaining a result of the 

correspondence application between images without 

resulting in obtaining any local solution when applying a 

correspondence between image. 

Furthermore, this invention provides a method capable 
25 of rapidly processing with a repetition process by each of 



blocks . 

Furthermore, this invention provides a method capable 
of performing a process considering over a distortion 
(deformation) of an input image. 
5 Furthermore, this invention provides a method capable 

of obtaining a similarity degree image being robust for 
variance of contrast of an image. 

Furthermore, this invention provides a method capable 
of widely using an information and reducing a case resulting 
10 in obtaining a local minimum. 

Furthermore, this invention provides a method capable 
of performing a process in which more widely information 
is used. 

Furthermore, this invention provides a method capable 
15 of performing an optimized minute adjustment by a repetition 
process . 

Furthermore, this invention provides a method capable 
of easily specifying corresponding points. 

Furthermore, this invention provides a method capable 
20 of saving useless accumulation-addition and efficiently 
correcting a deformation. 

Furthermore, this invention provides a computer 
program capable of realizing the method according to the 
present invention on a computer. 
25 Although the invention has been described with respect 



to a specific embodiment for a 
the appended claims are not 
be construed as embodying all 
constructions that may occur 
5 fairly fall within the basi 



complete and clear disclosure, 
to be thus limited but are to 
modifications and alternative 
to one skilled in the art which 
c teaching herein set forth. 
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